Bath generated work extraction and inversion-free gain in two-level systems 
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The spin-boson model, often used in NMR and ESR physics, quantum optics and spintronics, is 
considered in a solvable limit to model a spin one-half particle interacting with a bosonic thermal 
bath. By applying external pulses to a non-equilibrium initial state of the spin, work can be ex- 
tracted from the thermalized bath. It occurs on the timescale T2 inherent to transversal ('quantum') 
fluctuations. The work (partly) arises from heat given off by the surrounding bath, while the spin 
entropy remains constant during a pulse. This presents a violation of the Clausius inequality and 
the Thomson formulation of the second law (cycles cost work) for the two-level system. 

Starting from a fully disordered state, coherence can be induced by employing the bath. Due to 
this, a gain from a positive-temperature (inversion-free) two-level system is shown to be possible. 
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After E.L. Hahn discovered the spin-echo in NMR 
physics pj, it was soon suspected to violate the second 
law for a recent discussion discussion see e.g. [§,1). 
As a precursor to this question, we shall investigate in 
this Letter whether single or double pulses on single spins 
coupled to a bath already mark a violation of this law. 

Recently we analyzed the thermodynamics of the 
Caldeira-Leggett model for a quantum harmonic oscil- 
lator coupled to a quantum harmonic bath [||||. At low 
temperatures various formulations of the second law are 
violated: the Clausius inequality dQ < TdS is broken, 
the rates of energy dispersion and entropy production 
can be negative, and certain cycles are possible where 
heat extracted from the bath is fully converted into work 
("perpetuum mobile"). These findings are nevertheless 
in agreement with the Thomson formulation of the sec- 
ond law (cycles cost work) applied to an equilibrium ini- 
tial distribution, for which an exact proof exists B . 

The cause of the breakdown the universal thermody- 
namic picture is the occurrence of a cloud of interaction 
modes (photons or phonons) around the central system. 
Such a cloud is already familiar from the dressed electron 
picture in quantum electrodynamics. For systems of our 
interest, the cloud is not present at high temperatures, 
but, due to the non- vanishing coupling to the bath, it 
builds up at low T, inducing non-thcrmodynamic physics. 

Two level systems like electrons or two-state 

atoms in a field, or a two level Joscphson junction, dis- 
play relaxational behavior due to coupling to a bath. 
Transversal and longitudinal correlations relax on time- 
scales T 2 and T\ , respectively, with typically T 2 <C 71 [|| , 
since energy transfer is not involved in the T 2 - process. 

The Hamiltonian of the problem reads: 
H = Hs + Hb+Hi, Hs = -a z + —a x , (1) 



Hb = ^2,7kj k a\ah, Hi = - y^fffc(aj + Qfc)^- 

k k 

This is a spin i interacting with a bath of harmonic os- 
cillators (spin-boson model 0j9|); Hs, and Hi stand 
for the Hamiltonians of the spin, the bath and their inter- 
action, respectively. a x , o y and a z = —ia x a y are Pauli 
matrices, and a k and are the creation and annihilation 
operators of the bath oscillator with the index k 7 while 
the gk are the coupling constants. For an electron in a 
magnetic field B, e — gusB is the energy, with g the 
gyro-magnetic factor and fis the Bohr magneton. We 
shall restrict ourselves to the model with A = 0, which 
is a prototype of a variety of physical systems 0, and 
known to be exactly solvable [j^D , since z-component of 
the spin is conserved, and with it the spin energy. Phys- 
ically it means that we restrict ourselves to times much 
less than T\. In NMR (ESR) physics g] the model repre- 
sents a spin (electron) interacting with a bath of phonons. 
In quantum optics it is suitable for describing a two-level 
atom interacting with a photonic bath. 

Starting from general physical arguments , one typ- 
ically takes the quasi-Ohmic spectral density of the bath 

JH = E#^- w ) = - eWr ' ( 2 ) 

where g is the dimensionless damping constant and the 
exponential cuts off the coupling at w > T, the maxi- 
mal frequency of the bath. As usual, the thermodynamic 
limit for the bath has been taken here. 

Since A = 0, one has conservation of a z [t) = a z (0) (in 
the Heisenberg picture) . The dynamics of the annihila- 
tion operator of the fc'th bath mode reads 

&*(*) = e-^ l a k (0) + fj£* (e^ 4 - l) . (3) 
This implies 



1 



k 

where we denoted the quantum noise operator 

m =53fffc[4(°) e<w * t +a*( ) e " fa '* t ]. ( 5 ) 

which will act as a random force on the spin, and where 



(6) 



showing that 1/T is the relaxation time of the bath. 

Separated initial state. To describe situations, where 
the spin was suddenly brought into the contact with 
the bath, e.g. an electron injected into semiconduc- 
tor, atom injected into a cavity, or exciton created by 
external radiation, we make the assumption that ini- 
tially, at t = 0, the spin and the bath are in a sep- 
arated state, the latter being Gibbsian at temperature 
T = 1/(3: p{0) = p s (0) ® exp(-(3H B )/Z B , where p s (0) 
is the initial density matrix of the spin. In this situa- 
tion the quantum noise is stationary and Gaussian with 
average zero and time-ordered autocorrelation function: 
Kr{t-t') = (T[r)(t)fj(t')]), where T stands for the time- 
ordering operator and the brackets for the trace over the 
initial state. For t > it holds that 



K r {t) = K{t) -ihG{t) = h 2 [t{t) -iG x {t)] 
where an explicit calculation yields 



m - r(i + £-if)r(i + x + if) 



Gx(t) 



-Tt - 7(i), 1 (t) = - arctanrt 



7T TT 

The spin operators b± — b x ± i b y satisfy 

b ± = % - [±e + fj(t)-G{t)]a ± 
n 

and have, with loq = s/h, the solution 
<r ± (t) = exp (± iu/ t ) II± (t, 0) <7± (0) 



ft±(*i,*o) = e - lGl[tl - to) T exp 
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ds ry(s) 



to 



(7) 

(8) 
(9) 

(10) 

(11) 
(12) 



Depending only on afe(0) and at(0), it commutes with 
o"x,y,x(0). Thus one gets 



(5- + (t))=exp(twot) (n+(t,0) )<o-+(0)), 



(13) 



Evaluating the time ordered product term-by-term with 
help of Wick's theorem and resumming, we obtain 



The result is stationary, as it should. Substituting fll4| ) 
into dl3|), the and real and imaginary parts give 

(b x (t)) = [cosw t (b x (0)) - smcu t (b y (0))} e"^ (15) 
(a y (t)) = [sinwoi (b x (0)) + cos uj t (b y (0))} e~t® (16) 



Fort> 1/r Eq. (§) brings 



t 



To 



1 7i 

5 r 



(17) 



T2 can thus be identified with the transversal decay time. 
For small g this is a strong enhancement of the quantum- 
timescale Ti/T. The most optimistic case, T2 ~ 100 s at 
room temperature, involves g ~ 10~ 16 , truly small. 
The density matrix of the spin reads 

Ps = g 1 1 + (**(*)) + (*»(*)) *w + (**(*)) ] ■ (l g ) 

Its von Neumann entropy equals SVn = — trpslnps = 
— pi lnpi — P2 hrp2, where p^a = \ ±k|(ct)|. In the course 
of time |(<?(t))| decays to |(ct 2 (0))|, which makes the von 
Neumann entropy increase. Since there is no heat flow 
- the energy is conserved - this is in agreement with a 
formulation of the second law: the entropy of closed sys- 
tem, or of an open system without energy transfer (the 
spin in contact with the bath), cannot decrease. 

A sudden pulse. We mostly consider the Hamilto- 
nian ([[]) with A = 0. A fast rotation around the x- 
axis is described by taking A ^ during a short time 
61; this is called a fast pulse If A ~ is 

large, the evolution operator describing the pulse be- 
comes Ui = exp(— i8\H{A)/K) w exp(^i 9b x ) + 0(Si), 
where 6=—6iA/K is the rotation angle, 



U 1 1 & y U\ = b z sin 9 + b y cos 9, 
[/-f 1 b z U\ = b z cos 9 — a y sin 9. 



(19) 
(20) 



n+(ii,*o) ) = exp[-£(ti -t )], 



(14) 



When suddenly switching A on and off, the state of the 
system does not change, so p(t + Si) — U\ p(t) U-j - . The 
work done by the source is the change of the total energy 
which reads, since [Ui,Tl(A)} = 0, 

Wi(t) = tr [p(t)(H(A) -H) + P (t + 6 1 )(H- H(A))} 

= trp(t)(Ui 1 HU 1 -H), (21) 

The work done in this variation appears to be 

Wi= -|(1- cob fl)<5-,(0))- I sine (5- w (i)) (22) 

+ 1(1- cos 6)G{t)- T ^ sine (o- B (t)) 

Our main interest is work extraction from the bath. In 
order to ensure that the pulse does not change the energy 
of the spin, we first consider the case e = 0, where the 
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spin has no energy. For small g, 6 
one gets 



-tt/2 and t > 1/T 



2tt 



**(0)) e-*/ T2 



(23) 



If for a fixed t, temperature is neither too large nor 
too small, Te~*/ T2 > HT/n, work can be extracted 
(Wi < 0), provided the spin started in a coherent state 
(a x (0)) = —1. This possibility to extract work from the 
bath disappears on the timescale T 2 , because then the spin 
looses its coherence, (a XiV (t)) — > 0. Without a pulse the 
spin energy is conserved, however. Notice that any com- 
bination of ±7r pulses (this is a classical variation, since 
the coherence is not involved) can extract work only from 
a non-thermalized bath, i.e. for times ~ 1/r. 

Initial preparation via a rotation. Our approach 
also allows to consider a specific, well controllable non- 
equilibrium initial state: a Gibbsian of the total system, 
Pg = exp(— pTl)/Z, in which at t = the spin is rotated 
over an angle — around the y-axis, p(0) = Uq pc U^ 1 , 
with Uq = exp(— iira y /4). This maps a x — > a z , a z — » 
—<j x . Such a state models the optical excitation of the 
spin, as it is done in NMR and spintronics. Though p(0) 
does not have the product form, the problem remains 
exactly solvable. Taking 9 = — one now gets 



W l 



G(t) 



£ , . . Be . 

— (sin luqI cos 7 tanh — + cos w i sin 7) 



-I — — (cos wq£ cos 7 tanh — sm wq£ sin 7) 



ghT 
2tt 



— sm woi H — — cos wq£ 



tanh^ e- 1 '^ 
2 



(24) 



where j(t) of Eq. (^|) arises from friction, with 7(00) = 
^g, and where wo = e/h. Typically g is small, so work is 
extracted (W\ < 0) when the sinus is positive. The work 
decomposes, 

W\ = AU — AQ (25) 
into the change in spin energy due to the pulse, 
AU= £ -[{a z {t + ))-{a z {t-))] = £ -{a x {t-)) 



- — [sin uj t cos 7 tanh — + cos ui t sin 7] e 

£ . , /3e _ f/r 

-— sin wn£ tanh — e ' 2 , 
2 2 



(26) 



and the heat absorbed from the bath 



hTcosw t tanh^-e~ t/j!2 

7T 2 



(27) 



Notice its similarity with —W\ of Eq. (|23|). An inter- 
esting case is where work is performed by the total sys- 
tem (Wi < 0) solely due to heat taken from the bath 



(AQ > 0, AU = 0). This process, possible by choosing 
t w 27rn/wo with integer n, can be considered as a cycle 
of a perpetuum mobile, forbidden by folklore minded for- 
mulations of the second law. Indeed, under a rotation the 
length I (a) \ , and with it the von Neumann entropy, is left 
invariant, so one has a process with AQ > , ASVn — 0; 
which violates the Clausius inequality AQ < TAS v n- 
The work needed at time zero to rotate the spin is 



Wo 



-trp G - 



Va 



e , Be qfiT , . 

2 tanh T + ^T (28) 



representing the work done on the spin and on the bath, 
respectively. It can be verified that the total work 
Wq + Wi is always positive, so Thomson's formulation for 
a cyclic change H (here: the combination of the pulses 
at time t = and t) starting from equilibrium is obeyed. 

Two pulses in a rotated initial Gibbsian state. If 
there are many spins, each in a slightly different external 
field, there appears an inhomogeneous broadening of the 
u>0 = e/h line, for which we assume the distribution 



p(^o) 



7T (wq ~ W ) 2 + [T 2 * 



(29) 



having average u)q and inverse width , typically much 
smaller than T 2 . In this case the gain for a single pulse is 
washed out, leaving only the loss AQ = —ghT /2n : so two 
pulses are needed. We consider again the rotated initial 
Gibbsian state, and perform a first — \tt pulse around the 
x-axis at time t\ and a second i-7r pulse at time t 2 = ti+T. 
In the regime of small g and large t\ ~ T 2 the work in 
the second pulse is 



w 2 = 9.K _ I e -vr 2£si tanh ^£ 

2tt 2 ° 2 



(30) 



1 Be 
——e~ t2 ^ 7 ' 2 tanh — coswo^i(e sinc^r + <?Tcos loqt) 

At moderate times only slowly oscillating terms survive. 
They are the ones that involve At = t 2 — 2t\. For the 
total work W\ -I- W2 this brings 

W= g —~ ^/r^-m/r; tanh ^0 { A 

7T 4 2 1 

1 sg(At) /3S,w _ 

(1 + . , )j cos w A^ } (31) 



1 9 



sinh BhujQ ' 



For At near 27rn/wo such that the odd terms cancel, this 
again exhibits work extracted solely from the bath. 

Bath-induced gain without inversion. It is common 
knowledge that a two-level system with population in- 
version, i.e. with a negative temperature, is capable to 
amplify light, and it represents the basic working mecha- 
nism of lasers and masers (see for more non-standard 
mechanisms of lasing). In this context a bath is typically 
considered to be a drawback for the amplification, as be- 
ing a source of undesirable noises and relaxation towards 



3 



equilibrium |L(|. Our present aim is to show that the 
bath can nevertheless play a totally different role, namely 
in assisting work extraction (gain) by means of an positive 
temperature state in the two-level system. This effect is 
strictly prohibited by the second law applied to a positive 
temperature spin state if there were no coupling to the 
bath Q. We consider separated initial conditions with 
(<Tx(t)) — (&y(t)) = 0; an d apply a — \~k pulse at time 
t\ = + , and a \tt pulse at t 2 . For t 2 ^> 1/r the work is: 



W = AU-AQ, 

AU = -~ [1 - e -« ta > cos^ t 2 ](«T z ) + ^ e-« t2 



Ag 



— + ^gTe-^ smu t 2 {a z ) 

7T Z 



(32) 
sin LOot 2 

(33) 



where £ was defined in (|l7|). In the inversion- free case, 
the initial state of the spin is a Gibbsian connected to 
a positive temperature T = f° r which (a z ) — 

— tanh ^/3 e < 0. Let us first investigate the case T — oo 
(completely random state, (a x ,y,z) — 0). The work W 
can be negative (gain) provided e > AhT/ir. This situ- 
ation can be met in quantum optical two-level systems 
0,0 and in NMR Q. This mechanism concerns work 
extraction with help of the bath (it disappears for g — > 0) , 
but not from the bath, since now AQ < 0. The origin of 
the effect is that although the state of the spin was com- 
pletely disordered initially, the first pulse does generate 
some coherence. Due to back- reaction of the bath one has 
after the pulses (a v (t 2 )) = sin7(i2) exp(— £(£2)) sinwo*2, 
where 7(^2) of Eq. (^) goes from to \g on the timescale 
1/r, the reaction time of the bath. 

At finite To the term AU can still be negative when 
Tq > e/g, which can be met for not-too-small g, a condi- 
tion anyhow needed for having a sizeable effect. ^From a 
thermodynamic point of view the gain can be seen as a 
flow of energy from a high temperature (of the spin) to 
a lower one (of the bath) , and the outside world (gain) . 

Feasibility. Let us present several reasons favoring the 
feasibility of the proposed setups: 1) Two-level systems 
are widespread, not the least because many quantum sys- 
tem act as two-level under proper conditions; 2) Detec- 
tion in these systems is relatively easy, since already one- 
time quantities (&(t)) completely determine the state; 3) 
The harmonic oscillator bath is universal 4) Work 
and heat were measured in NMR experiments more than 
35 years ago ; 5) Our main effects do survive the av- 
eraging over disordered ensembles of spins, thus allowing 
many-spin measurements. 6) The ongoing activity for 
implementation of quantum computers provides experi- 
mentally realized examples of two-level systems, which 
have sufficiently long T 2 times, and admit external vari- 
ations on times smaller than T 2 : (i) for atoms in opti- 
cal traps T 2 ~ Is, 1/F ~ 10~ 8 s, and there are efficient 
methods for creating non-equilibrium initial states and 
manipulating atoms by external laser pulses H]; (ii) for 
an electronic spin injected or optically excited in a semi- 



conductor T 2 ~ 1 /us |15[ ; {Hi) for an exciton created in 
a quantum dot T 2 ~ 10~ 9 s Q (in cases (ii) and (Hi) 
1/r ~ 10~ 13 s and femtosecond (10 _15 s) laser pulses are 
available); (iv) in NMR physics T 2 ~ 1 — 10s and the 
duration of pulses can be comparable with 1/F ~ 1 /is. 

Conclusion. Like in the oscillator model |I| , the Clau- 
sius inequality can be violated in the spin-boson model, 
namely by a single pulse. For many spins this typically 
does not happen in the first pulse; it may occur in the 
second pulse, while it is absent for the two pulses taken 
together, because of intermediate entropy increase. Work 
can be extracted from the equilibrated bath, for the sin- 
gle spin case by one pulse, and for the many spin case by 
two pulses. This contradicts Thomson's formulation of 
the second law: no gain from cycles. Gain is also possible 
from a positive temperature (inversion-free) initial state, 
which may serve as a principle for bath generated lasing 
and masing. These effects are not in a conflict with the 
equilibrium Thomson formulation of the second law, and 
are fully quantum: they disappear at high temperatures 
and they do not occur for classical pulses. 
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